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Abstract 

The functional Schrodinger equation (FSE) for QCD gives a unique 
perspective on generation of a gluon mass m, as required for center 
vortices. The FSE, which yields a special d=3 gauge action, com- 
bined with lattice calculations strictly in d=3 give a value for the 
dimensionless ratio of d=3 coupling to mass g\jm. From this we in- 
fer a reasonably accurate value for the d=4 running coupling <? 2 (0) in 
the region of low momentum where it is nearly constant. The result, 
consistent with other estimates, is too low to drive chiral symmetry 
breaking (CSB) for quarks in a standard gap equation that has no 
explicit confinement effects. We recall and improve on old work show- 
ing that confinement implies CSB for quarks, and consider CSB for 
test (that is, quenched) Dirac fermions in the adjoint representation. 
Here the previously- found value of g 2 (0) is large enough to drive CSB 
in a gap equation, which we relate to the presence of center vortices 
(non-confining, for the adjoint) and nexuses that drive fermionic zero 
modes. We discuss the extension of adjoint CSB to finite temperature. 
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1 Introduction 



I discuss two related topics. The first is an approximate evaluation of the 
QCD coupling at zero momentum, a s (0), using the functional Schrodinger 
equation (FSE). The second is the beginning of a program, not completed, for 
relating chiral symmetry breakdown (CSB) properties for both quarks and 
(hypothetical) adjoint Dirac fermions in QCD to a s (0). This program can 
be looked on as using known CSB results (from the lattice) to constrain the 
allowed range of a s (0) or as using a theoretically-determined value of a s (0) 
to predict the CSB properties. The present accuracy of some preliminary 
investigations of these issues is not high enough to be definitive, but it does 
suggest that a s (0) in the range 0.6±0.2 is consistent with known properties 
of CSB. 

The FSE has long been of interest to me pQ; more recently, motivated by 
a talk [2] of S. Olem'k at the 2006 Oberwolz symposium, I wrote a paper last 
year [3], one conclusion of which is an estimate of the d = 4 strong coupling 
constant^ at zero (or small) momentum, a s (0) = g 2 (0)/(47r) ~ 0.5. This 
estimate is somewhat unusual since this estimate is based on one-loop gluon 
gap equations and lattice numerics in d — 3, not d = 4, plus some theory 
that attempts to relate d = 3 QCD to d — 4 QCD through the FSE. The 
d = 3 studies [niinilllllSlIISlinillSlIiaEniEI] vary somewhat, but in my 
view (Ref. [16J argues that the one-loop gap-equation results for m are too 
small) the most reliable value for the ratio g\jm is ~ 6.3/iV (for gauge group 
SU(N)), where m is the dynamical gluon mass and g\ is the d = 3 gauge 
coupling. FSE theoretical estimates suggest g 2 (0)/(47r) ~ (2.9/47r)(g 2 /m) ~ 
0.5. There are, of course, other estimates of az s (0): Phenomenology sensitive 
to infrared properties of QCD gives a s (0) ~ 0.7 ± 0.3 [22J. Other pinch- 
technique calculations [9], [10] suggest a s (0) ~ 0.5, just as I use here. 

Recently, lattice simulations [23] have been reported for an interesting 
question with a long lineage [2U, [25J [26]: How does CSB work in QCD with 
adjoint fermions^) instead of the usual fundamental-representation quarks? 

2 Please note that this coupling is the scheme- and process-independent coupling defined 
in the pinch technique [U [10] and not the process-dependent running coupling 
that is related to our coupling by a process-dependent transformation [11]. The pinch 
technique is an all-order way of extracting from Feynmans graphs for the S-matrix new 
graphical structures for off-shell proper vertices that are completely locally gauge-invariant. 

3 The lattice work is, of course, in Euclidean space, where there are no Majorana 
fermions and hence adjoint fermions plus gauge fields are not a supersymmetric theory. 
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These authors' works state that there is CSB for quenched adjoint fermions 
not only at temperatures below the deconfinement phase transition but also 
up to a temperature exceeding the deconfinement phase transition. In con- 
trast, for quarks there is no CSB above the deconfinement phase transition, 
consistent with lattice findings [27] that center vortices, the standard con- 
finement machinery of today, are both necessary and sufficient for CSB with 
quarks. Adjoint fermions, blind to the long-range parts of center vortices, are 
not confined and the argument is, therefore, that there are different mecha- 
nisms for CSB for fundamental and adjoint fermions. I do not claim to have 
definitively answered this question, but present numerical estimates based 
both on estimates of a s (0) and on fermion gap equations that have appro- 
priate kernels for small momentum suggest an answer that, given fairly large 
uncertainties, seems to accord with present-day knowledge from the lattice. 
What I report here is at best the beginning of a program of refining our 
quantitative understanding of QCD by comparing estimates of a s (0) with 
lattice and theoretical studies of adjoint CSB. 

Theoretical papers from long ago ESI ED] argue that various con- 
finement mechanisms produce CSB for quarks. (Later I will give a brief 
update of some of these ideas, based on the role of a condensate of cen- 
ter vortices and their close relatives, nexuses, in confinement and CSB). 
These confining mechanisms depend only indirectly on a s (0). In princi- 
ple it could be (although the lattice data say otherwise) that there would 
be CSB for quarks even with no confinement, a possibility that does de- 
pend on a s (0). Generally, CSB for quarks and adjoint fermions should be 
sensitive to three couplings: (1) the standard QCD coupling g 2 (0)/(47r) at 
zero momentum transfer; (2) the critical coupling g 2 (0; fund) /(47r) above 
which CSB occurs for fundamental-representation fermions (quarks) as found 
from a gap equation that does not contain confinement effects; and (3) 
g^(0; adj)/ (4n), the same coupling for (quenched) adjoint fermions. Because 
these last two couplings differ from the first only by Casimir factors, the criti- 
cal couplings for gap equations are inversely proportional to C2, the quadratic 
Casimir eigenvalue for the fermions in the gap equation, and so for QCD, 
<7c(0; fund) / (An) = (9/4)^(0; adj)/(4<K). However, just knowing this is not 
enough to settle the issue of whether CSB can or cannot take place through 
a standard gap equation for quarks or adjoint fermions; we need not just the 
ratios but also the values of the critical couplings to compare with estimates 
from other QCD models of the couplings. 
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There are in principle three possibilities: 



9 2 (0) 
g*(Q;adj) 
g 2 c (0; fund) 



< 



< 



< 



gl(0;adj) 

g 2 (0)<g 2 c (0-Jund) 
9 2 (0) 



(1) 



In the first case only quarks can show CSB, and confinement is necessary 
for quark CSB. In the second case (which is favored both from the lattice 
data and from the estimates I give here) there is CSB for both adjoint and 
fundamental fermions. Confinement is again necessary for quark CSB, whose 
transition temperature must be rather near that of the deconfinement transi- 
tion; the gap equation is largely irrelevant, althouogh it may account for some 
separation in these two transition temperatures.. CSB for adjoint fermions 
comes solely from non-confining effects as summarized in a gap equation, and 
adjoint-fermion CSB may or may not extend above the deconfinement tran- 
sition. In the third case, there is CSB for both kinds of fermions, and this 
ought to persist even if confiners such as center vortices are removed from 
the lattice simulations. Since there is more than a factor of two between 
adjoint and fundamental critical couplings, the inequality in the second case 
is a fairly broad one and great accuracy is not needed to single out this case. 

Fermion gap equations have a long history, beginning with the JBW equa- 
tion [21], and nearly all previous work that does not address confinement 
issues makes three approximations: 1) Gluons are massless; 2) Landau gauge 
is used; 3) because to one-loop order in this gauge vertex corrections are 
not ultraviolet-divergent, vertex corrections are ignored. Non-perturbative 
phenomena of low-energy QCD require a more careful treatment that I will 
sketch here. In the first place, infrared slavery implies dynamic gluon mass 
generation, which can be only studied effectively in the gauge-invariant pinch 
technique [61 [9] . In the second place, because fermion mass generation is an 
infrared effect there may be important low-energy fermion-gluon vertex cor- 
rections. Finally, the possibility of such corrections requires a more careful 
study of gauge invariance of the gap equation. To some extent these issues 
have been addressed before [321 EH Ell > but not in a context particularly use- 
ful here. In the present work I include a dynamical gluon mass in the gluon 
propagator; give a sketch of the derivation of a gauge-invariant gap equation 
using the pinch technique (which seems not to have been addressed in detail 
before); and use the gauge technique to infer approximate low-momentum 
vertex corrections that satisfy the correct Ward identities of the pinch tech- 
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nique. The simplest application of these principles, all that I report here, 
yields a gap equation much like the JBW [3T] equation except that the gluon 
is massive. 

The work reported here is still very much in progress, and needs consid- 
erable sharpening. However, I believe it is already at the point where one 
can qualitatively see and explain the differences between CSB for quarks and 
for adjoint fermions. 

2 The functional Schrodinger equation for QCD 

This work has already been published [3] so I will be brief here. The vacuum 
wave functional of QCD is a gauge-invariant functional that I write in the 
form: 

V>{A a (z)} = e- S:i{A ^ £)} (2) 

in which 2S% is a real gauge-invariant d = 3 effective action (a factor of 
two because \ip\ 2 is the weight function for constructing vacuum expectation 
values). It is constructed to satisfy 

Hif> = E vac ij (3) 

with Hamiltonian 

H = J {-^(^) 2 + ^(^) 2 ) - / \\ W) 2 ] + V. (4) 
The functional S3 has infinitely many terms: 

g 2 Ss = \\ J A'Hkr\" + + • • • (5) 

Gauge invariance requires [TJ that Qij be conserved, and any two successive 
terms in the expansion are related by ghost-free Ward identities. Even the 
simplest of these terms for a free gauge theory is not familiar as an effective 
action, because it has a square root: 

Ssfree = ^ / A?V=V*P^. (6) 

Here Pij is the usual transverse projector: 



My approach to the FSE (and later to the fermion gap equations) begins 
with the fact that QCD, because of infrared slavery, undergoes dynamical 
gluon mass generation jH (6J [7J [131 HO]. The FSE must display this fact; 
how does it do so? It is easy to describe in a toy model, an Abelian gauge 
theory with gauge-invariant mass M put in by hand. The Hamiltonian is: 




The corresponding S3 that exactly satisfies this Hamiltonian is: 

S 3A bei = ^2 / ^Wm 2 - ^PijAj. (9) 

Once I add the mass, a nonlocality appears from the transverse projector. 
But this is easily remedied, by introducing a scalar field in the mass termQ 

in f 

S 3m Abel = ^2jl A i- ^ ( 10 ) 

and functionally integrating over not only the gauge potentials but also over 
when constructing vacuum expectation values. This is entirely equivalent 
to using the non-local transverse projector. Note that the mass term by itself 
satisfies a FSE with the Fj* term missing from the Hamiltonian of Eq. (jSJ). 

Since we are interested in infrared effects, it is reasonable to make an 
expansion in inverse powers of the mass m (or equivalently in powers of the 
gradient operator). The leading term is 0{m) and is just the mass term 
SzmAbei itself. However, a naive expansion runs into a little bit of trouble, as 
we see by expanding the square root in the exact Abelian solution: 

S.Abel =-2ji/ ~ 2^ + • • ' ]PijAj - (U) 

(Observe that the second and succeeding terms in the expansion are local; in 
fact, the second term is, up to a factor 1/m, the usual Ff- term.) If only these 
two terms are kept, the field has mass \f2m and not m. Higher-order terms 
not written must correct for this discrepancy. Of course, no such expansion 

4 Or one can begin with a simple gauge-dependent mass term J mA 2 and project out 
its gauge-independent part by integrating over all gauge transformations of the gauge 
potential. 
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in 1/m is necessary for the Abelian case, but it is for the non-Abelian case, 
and requires [3j a somewhat better approximation to the square root which 
is usable for momenta whose components are comparable in size to m. 

How does this go for the non-Abelian case? I have argued for decades 
[35] that one describes locally gauge-invariant gauge-boson mass generation 
through a gauged non- linear sigma model (GNSM), analogous to the local 
mass action S^m above. To simplify the notation I use the anti-Hermitean 
matrix gauge potential 

A = ±-.KAt, (12) 

where the A a are the standard Gell-Mann matrices, and covariant derivative 

Di = di + Ai. (13) 

Introduce a unitary matrix U, with the gauge transformation properties U — > 
VU when the gauge potential is transformed by: 

Ai -> VAiV- 1 + Vd l V-\ (14) 

Then the locally gauge-invariant GNSM mass term i^j: 

S 3m = Z ^ [ Trp-'D.U} 2 . (15) 
9 J 

The non-covariant derivative UdiU' 1 is the non-Abelian generalization of the 
Abelian scalar <9,0; in fact, the GNSM action can be written as: 

S, m = ^ [ Tr{A - UdtU- 1 ] 2 . (16) 

One can, just as in the Abelian case, eliminate U through its equations of 
motion (that is, minimize S , 3 OT ), which are (after some non-trivial algebra): 

[D i ,A i -Ud i U- 1 ]=0. (17) 

The perturbative solution has infinitely many terms, of which a few are [35J: 

U = e^ u = ^ d .A+±{[A i A^d.A]+±[d.A,±d.A) + ...} 

(18) 

5 Since I am only interested in infrared effects I intepret the coupling as being evaluated 
at zero momentum. 
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The linear term simply generates the transverse projector I have already 
used in the Abelian case. In addition, there are non-perturbative solutions 
relevant for center vortices. 

I claim that this GNSM mass term is the leading term in the 1/m expan- 
sion of an effective action that capture the leading non-perturbative effect of 
infrared slavery, which is dynamic gluon mass generation. (It also captures 
the structure of massless poles in the pinch-technique Schwinger-Dyson equa- 
tion yielding the mass dynamically.) It is a good candidate for the leading 
term of the effective action S3 because it is gauge-invariant; indeed, just as in 
the Abelian case, it comes from projecting out the non- Abelian gauge invari- 
ant from the simple A 2 mass term. It is almost evident without calculation^! 
that the next-leading term should be the usual Yang-Mills term, which is just 
the non- Abelian gauge completion of the Abelian term shown in Eq. (fTTj) : 

S 3 = ^ / Tr[Ai - UdtU- 1 } 2 - -L- ( TrG% + ... (19) 

where G^- is the usual Yang-Mills field strength. 

The normalization of the second term follows from the fact that the 
quadratic term in S3 is just the Abelian action of Eq. pip , one copy for 
each gauge boson. But just as in the Abelian case this wrongly yields a 
free-field mass of \/2m instead of m. I have proposed (3J to cure this ap- 
proximately by choosing a renormalization factor Z that best approximates 
the actual square root operator of the FSE with a two-term expansion; the 
result for the approximate two-term action 7^=3 = 2S3 is: 

- 2S 3 = -7^3 = ^ / d'xTrp-'D.U] 2 } + SxTrG% + 0(m" 3 ) 

(20) 

with Z ~ 1.1- 1.2. 

The final step is to compare this to the standard d = 3 form of gauge- 
invariant massive QCD for the given mass m, which is: 

I d=3 = -Jd 3 x \^TrG% + ^Tr^At/] 2 } • (21) 

Here the d = 3 coupling g\ has the dimensions of mass. Comparing the two 

6 Ref. j3] contains details about the calculation, using the pinch technique and the 
gauge technique. 
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forms of Ia=3 yields: 

9 2 = (22) 
m 

Several authors [T2llT3linilI5llISlinilIHllISll2Dll2I] have given either lattice 
or theoretical estimates of the dimensionless ratio mjNg\ for gauge group 
SU (N) with N = 2,3. The results are quite consistent with an average value 
of gl/m ~ 6.3/ N. One can roughly convert the no-quark coupling in Eq. (1221) 
to three light flavors by multiplying the right side of Eq. (122]) by 11/9, and 
the resulting value of the strong coupling at zero momentum is: 

^ * 0.5. (23) 

This estimate and other quite similar pinch-technique estimates [9j \10\ , 
combined with phenomenology [22] that gives somewhat higher values, sug- 
gests that a s (0) ~ 0.6 ± 0.2. I will now see how this range of values fits into 
fermion gap equations for CSB. 



3 Confinement, soliton condensates, and gap 
equations 

The first step is to understand the difference between gap equations that 
purport to show the effects of confinement and those that do not. I will 
not do that in any detail here, but simply draw a few conclusions from the 
fact that confinement comes from the long-range pure-gauge parts of center 
vortices, which are quantum solitons of an effective action of the type given in 
Eq. (}2~Tj) or its d = 4 extension. Because of complications having to do with 
integrating over center-vortex collective coordinates, it is easiest to present 
the argument in d — 2, where I will use not the familiar action of d = 2 Yang- 
Mills theory but rather the effective action, with a mass term, of Eq. (1211) in 
two dimensions [36J. (I have argued [3] that this is (possibly up to an overall 
factor) the correct action for the d — 2 + 1 FSE vacuum wave functional, and 
that leaving out the mass term cannot be right; the reason is that without 
the mass term Wilson loops of all representations show an area law, while 
with it only iV-ality ^ representations are confined, which is correct for 
d = 2 + 1. But all this is irrelevant to the present argument.) 
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The simplest center vortex is a soliton solution to the equations of motion 
of the effective d = 2 action: 

Aj(x -a;K) = {2v:Q K /i)e jk d k {A m (x - a) - A (x - a)}. (24) 

Here Qk, K — 1 . . . N — lisa generator of the Cartan subalgebra of SU (N), 
normalized so that exp[27uQi\;] is an element of the center, and A mi o are 
free propagators of mass m, 0. The vector a is a collective coordinate for 
translations, and I do not indicate collective coordinates for group rotations. 
A gluon propagator can be defined by integrating the product of two soliton 
potentials over their common collective coordinates: 

x W2] 

mx)A){y)) = -J*- £ £ (-2)Tr(Q 2 J )A i (x - a; K)A 3 (y - a; K). (25) 

N 1 a K=l 

This propagator has a long-range part coming from the Ao term, and 
the remainder is short-range. The full gluon propagator, with both terms, 
can be used in the gap equation for the pinch-technique fermion propagator. 
This gap equation is derived from an S-matrix element, or equivalently from 
some complicated functional of a Wilson loop. The long-range pure-gauge 
parts are detected by their linkage with the Wilson loop. If a vortex is inside 
the loop it gives a non-trivial center element; otherwise it gives unity. Since 
there are no non-trivial elements for iV-ality zero representations, such as the 
adjoint, the adjoint fermion is completely blind to them and sees only the 
short-range parts. So for quarks the gap equation should be described with 
a gluon propagator containing the long-range pure-gauge terms, which gives 
the long-range gluon propagator: 

(A^(x)A b j {y))\i ong = consUabiAa^x - y) (26) 

where 

(Ao)y(x -y) = -^y 2 J d 2 k{5 tJ - k^/k 2 )— (27) 

is the gauge propagator of d — 2 QCD. This propagator, singular at large 
distances, not only confines quarks, it breaks CSB. (Generally in d dimensions 
the propagator behaves like k~ d , which I used in d = 4 in an earlier discussion 
[28] of confinement and CSB.) 

The remaining short-range part couples to all fermions with strength 
proportional to the quadratic Casimir C*2, and has range 1/m. In fact, these 
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short-range soliton parts must sum to a standard massive gauge propagator 
of the form (omitting the group labels): 



Aijik) = Pij(k) 



k 2 + m 



1 




C,k{kj 



(28) 



I have written this propagator in the form given by the pinch technique, where 
the physical part is both gauge-invariant (unlike conventional propagators) 
and transverse. The last term on the right is a necessary but inert term de- 
pending on the chosen gauge that cannot enter any pinch-technique physical 
prediction. In particular, it must cancel out in pinch-technique fermion gap 
equations. 

4 Fermion gap equations without confinement 

In this section I briefly mention the older gap equations, which are oriented 
toward ultraviolet behavior. Then I go onto newer equations that treat the 
infrared regime of QCD more accurately, including a quick discussion of a 
pinch-technique gap equation. More details on these newer gap equations 
will be published elsewhere. 

The history of gap equations, from the Johnson-Baker- Willey (JBW) 
equation of the sixties [31] to work of the nineties, can be traced from various 
specializations of an approximate gap equation for the CSB-breaking running 
fermion mass M(p 2 ): 



An often-studied variant, and the only one I will consider explicitly, drops 
the non-linear fermion mass terms in the denominator of the fermion propa- 
gator on the right side of the equation, replacing k 2 + M(k 2 ) 2 by k 2 . Trouble 
arises with any linearized gap equation that has a massless gluon propagator, 
because it is impossible to have a finite fermion mass M(0) at zero momen- 
tum. Removing this problem by keeping the full non-linear equation with a 
massless gluon and a fermion mass in the denominator is not really a good 
solution, since the only mass scale would have to come from the running 
charge, which is essentially constant in the low-momentum regime where the 
non-linear fermion mass is needed as an infrared cutoff. In fact the funda- 
mental infrared cutoff is the gluon mass m and I expect that M(0) ~ m. 



M{p 2 ) = 3C 2 



I 



d 4 k g 2 (k 2 )M(k 2 ) 

(2vr) 4 (k 2 + M(k 2 ) 2 ){{p - k) 2 + m 2 ) ' 



(29) 
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JBW used the linearized equation for QED, where Ci = 1, the charge does 
not run, the photon mass m vanishes, and M is the electron mass, supposed 
to be generated spontaneously in massless QED. The idea was to show that 
this equation leads to a running mass vanishing at large momentum, hence 
requiring no bare mass counterterm. When this equation is used for QCD, 
C*2 is the quadratic Casimir eigenvalue, m is the gluon mass, and g 2 (k 2 ) is 
the running QCD charge. (In principle the dynamical gluon mass m 2 must 
run too [6], but since the most important effect of this mass is in the infrared 
region I will not include such running.) 

The master gap equation is a simplified form of the Schwinger-Dyson 
equations for the fermion propagator S(p), which has the form: 



As defined, the fermion "mass" is essentially gauge-invariant, at least in the 
sense that its ultraviolet anomalous dimension is gauge-invariant. In practice, 
most workers specialize to the Landau gauge because large-momentum radia- 
tive corrections to the fermion- gluon vertex are absent in one- loop order, so 
it is argued that in this gauge it should be a decent approximation to ignore 
vertex corrections and set (using the QED Ward identity) A(p 2 ) = 1 for all 
momentaQ I will also use Landau gauge and assume A(p 2 ) does not change 
much with momentum, although this requires a few words of justification for 
the pinch technique/gauge technique approach given later. 

The linearized fermion gap equation has the generic matrix form 



where K is the kernel, derived from the single skeleton graph for the inverse 
fermion propagator. There is always a chiral-symmetry-preserving solution 
M = 0, but we seek CSB-breaking solutions with M/0. If the kernel is a 
well-behaved (finite-dimensional, bounded) matrix it is clear that CSB can 
only occur if g 2 is sufficiently large; otherwise the determinant det(l — g 2 K) 
will not vanish. Actually, K (from, for example, Eq. (|29j) ) is not that well- 
behaved, but in the equations we use there is a critical coupling g 2 marking 
the boundary between CSB and chiral symmetry preservation; most students 
of gap equations give rather similar values for this critical coupling. 

'In QED, with its small coupling, ignoring higher-order effects could well be justified. 
In QCD, with its strong coupling, the justification is that we are looking for infrared- 
dominated effects, so large-momentum contributions should not be important. 



S-\p) = p ! A{p 2 )[l + iM{p 2 )}. 



(30) 



M = g 2 KM 



(31) 
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There are two forms of criticality: The first is based on a differential 
equation derived from the gap equation. Its solution may be well-behaved 
for sufficiently small coupling, and then show unphysical features, such as 
alternation of signs of the running fermion mass, for larger coupling. The 
second is based on the original linearized integral equation, which imposes 
a boundary condition equivalent to consistency between the left-hand and 
right-hand sides of the gap equation evaluated at zero momentum. Consis- 
tency fails if the coupling is too small. Whether either or both criticality 
criteria hold depends on the gap equation, as I will show by an explicit ex- 
ample. 



4.1 Ultraviolet behavior: The JBW equation and vari- 
ants 

Now I give some simple special cases, the first of which is well-known, of 
the master gap equation of Eq. (129]) . The first is the original JBW equation 
used to study possible dynamical generation of the electron mass in QED. 
The JBW equation sets g 2 to a constant and has no mass terms in the 
propagators on the right side. With the aid of 

nj -± w = -^S( P -k) (32) 

the JBW equation (now for QCD, so C*2 is reinstated) becomes the differential 
equation 

„ 2M' AM , , 

M H — H — = 33 

where 

and the primes indicate derivatives with respect to p 2 . There are two linearly- 
independent solutions: 

M ± (p 2 ) = const. (p 2 ) v± , u± = ± [1 - AX} 1 ' 2 } (35) 

For small coupling the V- solution decreases roughly at large momentum like 
1/p 2 and is commonly called infrared-dominated; the v + solution corresponds 
to the original JBW solution, which falls off very slowly at large momentum 
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and is called ultraviolet-dominated. Although we may term one of the solu- 
tions infrared-dominated, this does not mean that it can be used for small 
momenta; in fact, both M±(0) diverge. The infrared-dominated solution is 
only useful in the ultraviolet. 

These solutions M± are not appropriate for finding the ultraviolet be- 
havior in QCD, where the charge runs at large momentum. This issue was 
first clarified by Lane [37J, who showed that the master gap equation with 
a running charge and no masses exactly captures the ultraviolet behavior of 
the fermion running mass. At large momentum the gluon mass m can be 
dropped in the denominators, and the running coupling has its usual behavior 
(at leading order) 

^ * mkw> (36) 

where b is the leading term in the /3-function {(3 = —bg 3 + . . .) and A is the 
QCD mass. With the aid of the renormalization group Lane showed that the 
ultraviolet behavior appropriate to CSB is: 

2 (in P 2 r 3c 2 

M{p — ; a = ——. 37 

p z Ibn^b 

I will recover this behavior later in a JBW-like equation with both a gluon 
mass and a running charge. 

4.2 Critical couplings 

The simple JBW equation, with a non-running charge evaluated at zero 
momentum, has a critical coupling at A = 1/4, corresponding to a critical 
coupling g 1 of value 

2 

£=3CV (38) 

For couplings larger than critical the exponents v± become complex, with 
v*> = and the asymptotic solutions both decay and oscillate, for example, 
like ~ p^ 1 cos [In p]. There is certainly no reason to accept as physical a 
running fermion mass that alternates in sign. I have already noted that for 
a well-behaved kernel g 2 must exceed a certain value for CSB to take place, 
yet for the simple JBW equation criticality marks the onset of apparently 
unphysical behavior. So is criticality in the JBW equation at all related to 
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CSB? The general answer is yes, although CSB may require a somewhat 
different value from g 2 . 

There are in fact a few reasons to believe that g 2 c above is indeed close 
to the true critical coupling for CSB. For example, Miransky [38J and his 
collaborators [39] have studied the positronium Bethe-Salpeter equation and 
find that if g 2 /(Air) > 7r/4 tachyonic levels appear. They relate these to 
vacuum rearrangment and scale-breaking phenomena, of which CSB is an 
example, and compare the oscillatory behavior of equations like the JBW 
equation with supercritical coupling to the quantum-mechanical problem of 
fall into the center or to the behavior of solutions to the Dirac equation with 
supercritical QED charge Za > 1. In fact, the JBW differential equation is 
nothing but the radial Schrodinger equation at zero angular momentum and 
energy for an attractive potential V(p 2 = r): 

V(r) = ^. (39) 

This potential shows fall into the center if A > 1/4, which is just the critical 
coupling given above. The QCD critical coupling of Refs. [381 ES] would 
be 71/(402), not much different from the JBW critical value. Others (see 
[4"9] and references therein) claim that the JBW value is the critical coupling 
for the pion Bethe-Salpeter equation to admit a massless pion. So I will 
assume that a critical coupling deduced from the gap equation is close to, if 
somewhat below, the critical coupling above which there is CSB. 

5 Infrared gap equations, gauge invariance, 
and the pinch technique 

The renormalization group cannot say anything about the behavior of M(p 2 ) 
at low momentum, where (among other things) it becomes necessary to in- 
clude the effects of the gluon mass m not only on the propagator but also on 
the running coupling. I and others (see [3D], which has many references to 
other works) claim that in QCD there is a quasi- conf or mal infrared regime 
where the running charge g 2 {k 2 ) is only slowly changing with momentum. 
Long ago it was argued [6] that a decent approximation to the running charge 
at both low (Euclidean) and high momentum, with the right sort of two-gluon 
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threshhold, is: 



« (* > = Mn(t W) - <«» 

This quasi-conformal coupling runs very slowly at k -C m. Because higher- 
order terms are neglected this expression cannot be more than perhaps 10% 
accurate. In the ultraviolet region a wide range of values of A does not change 
the coupling very much, but in the infrared regime dimensional transmutation 
has taken place, with the zero-momentum coupling determined by the ratio 
2m/ A. The reader can verify that the choice m = 0.5GeV, A = 0.3 GeV 
is within 10% of a recent evaluation from data [41j of the strong coupling 
at the masses of the r and Z, and (for three light flavors) gives a s (0) = 
# 2 (0)/(4tt) ~ 0.6. 

Now consider the linear gap equation keeping the gluon mass, but the 
running charge is replaced by the fixed-point value g 2 {0). There is no simple 
differential equation, but one can do the angular integrals. The resulting 
integral equation is: 

M(p 2 ) = J™ d k 2 ^^M(k 2 )K(p, k) (41) 

where 

K(p, k) = K(k, P ) = p2 + k2+m2 + p + f£2+ m2 y _ 4p2fc 2]l/2 " ^ 

These equations can only be solved numerically, but they are closely re- 
lated to a solvable differential equation with a dominating kernel K, such 
that K < K. The new kernel is: 



The new kernel K is exactly equal to K for p > 0, k = (or k > 0,p = 0), 
and K, K are asymptotically the same at large momentum. When both 
momenta are non-zero, K is greater, by a maximum factor of about 1.3 at 
k 2 = p 2 = m 2 . It should therefore be that the critical value g 2 for K is greater 
than that for K. I expect a further (modest) increase in the critical coupling 
because if I had included a properly-running charge in the equation with 
kernel K it would further reduce this compared to the presently-considered 
kernel K with a fixed charge. 
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The approximate version of the original (linearized) master gap equation, 
using the kernel K and fixed charge, is: 

M(p 2 ) = 3 ^y\' f dk 2 M (k 2 ) + r J^-Mk 2 )} (44) 
yF ' 16tt 2 l p 2 + m 2 Jo K J h 2 k 2 + m 2 v ;/ v ; 

where (k, a positive number between zero and one, measures the discrepancy 
between using the kernel g 2 (k 2 )K and the kernel g 2 (0)K. One should think 
of (k as roughly measuring some sort of momentum average of the form 

^ K( P ,k)H4m 2 /A 2 ) 
X k{p,k)ln[{k 2 + 4m 2 )/A 2 } 

I estimate roughly that (k — 0.7 — 0.8. 

The differential equation emerging from Eq. (jUJ) is exactly the same as 
Eq. (1331) with no gluon mass, except that the independent variable is changed 
from p 2 to p 2 + m 2 and A of Eq. (151]) is changed to C&-A: 



2M' C K XM , . 

M" + ^— -+ ( I* 2 , 2 =0- 46 
p z + mr (p z + m A y 



The solutions are: 



M ± (p 2 ) = const.(p 2 + m 2 )^, u± = ^{-1 ± [1 - 4^A] 1 / 2 } (47) 

and, unlike the massless-gluon solutions of Eq. (1331) . these are finite at p 2 = 0. 
Criticality occurs now at Of A = l/4 ; or 

9 -M ^ (48) 

As mentioned above, there is another criticality criterion. I consider 
evaluating the integral in Eq. (T441) at zero momentum: 

1 CkX 



oo 

dk 2 K(k,0)M(k 2 )/M(0). (49) 

This sets a boundary condition on the linear combination of solutions M± 
of the differential equation, from which the solution of the integral equation 
must be formed. One might expect that if A is too small this criterion can 
never be satisfied, so the differential equation can be satisfied but not the 
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corresponding integral equation. However, because both M + and M_ are 
finite at the origin and integrable over the kernel K at infinity, it is in fact 
always possible to find a linear combination that satisfies Eq. f|49l) . The 
problem is to find two numbers a± such that 

a+M+(0) + a_M_(0) = 1 (50) 

and 

l = OcA r dk 2^^ k ^ a+ M+ (k 2 ) + a.M„(k 2 )}. (51) 
4 jo 

These 2x2 linear equations are soluble for a± except for at most one value 
of CkA. (In the case at hand, a± = v T /(v T — z/±); even the limit C^A = 1/4 
exists.) So the criterion of Eq. ( |49l) is not useful whenever the two solutions 
to the differential equation are both finite at the origin and integrable over 
the kernel at infinity. 

The final change to be made in the gap equation is to replace <? 2 (0) by 
g 2 (k 2 ) in the i^-equation given as Eq. (j4~4"l) . As long as the momentum 
dependence of the running charge is on the integration variable k 2 alone, a 
differential equation can be found for any fc-dependence. This differential 
equation is: 

2M' | 3C 2 C K g 2 ( P 2 )M 
p 2 + m 2 16ir 2 (p 2 + m 2 



M" + + yp ' = (52) 



Here (' K accounts for the average difference between K and K but not for 
the running charge (cf. Eq. f|45|) ). I estimate (' K ~ 0.9. Using the running 
charge of Eq. (|4"0"1) gives the equation: 

M" + 2M ' + - (53) 

p 2 + m? (p 2 + m 2 ) 2 ln[(p 2 + 4m 2 )/A 2 ] 



where a is Lane's constant, from Eq. ([3] 

To my knowledge this is not a form of any standard differential equation, 
but it reduces to one in two cases. For small momentum if p 2 is dropped 
compared to 4m 2 , the result is Eq. ( Ho*]) , already solved. An equation that 
is infrared-finite and asymptotically-exact for large momentum results from 
replacing p 2 + 4m 2 by A(p 2 + m 2 ), which yields a confluent hypergeometric 
equation. 

Unfortunately, this is not very accurate for small momentum, but since it 
is asymptotically-exact for large momentum I will consider it briefly. Replace 
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the running charge by: 



9 ^ k2) 61n[4(£; 2 + m 2 )/A 2 ]' (54) 

As one should insist, g 2 {0) is unchanged by this modification, and the large- 
momentum behavior is insensitive to m. 

In addition to the correction factor (' K , the last term of Eq. fl53l) should 
be multiplied by a factor ( g > 1 that attempts to correct for this mutilation 
of the running charge. Think of ( g as an average of the type: 

^ ln[4( P 2 + m 2 )/A 2 ] 
^ _ Mn[(p 2 + 4m 2 )/A 2j - 1 j 

Now change to the new independent variable 

With this factor, Eq. fl53|) becomes: 

M + M+ <9 ^ M =0 (57) 

where dots indicate i derivatives. A simple estimate suggests that ( g ~ 
1.1 — 1.2, so the product CgCx * s essentially unity within the accuracy to 
which I aspire, and I drop this product. 

Equation fl56|) is a confluent hypergeometric equation (see, for exam- 
ple, [12] )■ The two linearly- independent solutions corresponding to Whit- 
taker functions have different asymptotic behaviors. At large momentum the 
infrared-dominated solution goes like 

M«(P 2 )- ^[l + 0(^)] (58) 

which is just the behavior found by Lane [37]. The ultraviolet-dominated 
solution goes like 

M uv ^(lnp 2 )- a [l + 0(^)}. (59) 

Both of these solutions are finite at p 2 = and both are integrable over 
the kernel (including the running charge), so the same situation arises as 
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with the closely-related Eq. (1441) : Except for at most one value of a, there 
is always a linear combination of Whittaker functions that satisfies the zero- 
momentum consistency condition analogous to Eq. fT49~]) that holds when the 
charge does not run. 

Our other criterion for criticality is the onset of zeros in M(p 2 ). The 
confluent hypergeometric equation, Eq. floT|) . shows critical behavior at a = 
a c = 1 — not in the sense of singularities in the solution, but, as for the 
original JBW equation, if a > a c zeros of the mass function M(p 2 ) set in 

The criticality condition a c = 1 is not really dynamical; it depends only 
on the particle spectrum and gauge group (see Eq. (13"T|) ). In QCD with 
three light flavors I find for quarks that a = 4/9, but for (quenched) adjoint 
fermions, a = 1. It appears unlikely that the modified JBW equation could 
lead to CSB for quarks, but it might well do so for adjoint fermions, in view 
of all the approximations and uncertainties of our development. 

I can, for what it is worth, convert the condition a c = 1 into a criterion for 
the critical coupling, combining Eqs. (I37|40p by eliminating the /3-function 
coefficient b. The result is: 

= < ). 60) 

4tt 3C 2 l ln(4m 2 /A 2 ) ; y ' 

This estimate differs from that coming from the ultraviolet JBW equation by 
the factor in parentheses. This factor plausibly varies from 1 (at m ~ 3.7A) 
to 2 (at m ~ 1-4A), given uncertainties in both m and the effective value of 
A that works best for a one-loop approximation. 



5.1 The pinch technique/gauge technique gap equa- 
tion 

My final infrared gap equation is based on the gauge technique (Ref. [44] and 
references therein), in which one "solves" the Ward identity for the fermion- 
gluon vertex in terms of the fermion propagator. The gauge technique is 
combined with the pinch techniquejf) which is a way (already mentioned) of 
finding off-shell Schwinger-Dyson equations, such as the fermion gap equa- 
tion, that are locally gauge- invariant. These are not, of course, the usual 
Schwinger-Dyson equations; pinch-technique proper vertices and self-energies 

8 See [45] for a thorough discussion of the two-loop pinch-technique fermion self-energy. 
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Figure 1: The on-shell quark-gluon amplitude for finding the pinch-technique 
quark proper self-energy. The last graph (h) is an example of a new kind of 
graph emerging from pinching. 



contain contributions from parts of graphs naively unrelated to the proper 
vertex under study. As a result, the Ward identities of a non-Abelian gauge 
theory are modified; they are just the naive Ward identities of a QED-like 
theory. 

The pinch technique begins by setting up an on-shall S-matrix elements 
containing the off-shell Green's function of interest — see Fig. [TJ 

Parts of numerators of some of these graphs contain longitudinal momenta 
which, when they strike an elementary vertex, trigger an elementary Ward 
identity of the form 

k^^ = S-\k)-S-\p-k). (61) 
When an inverse propagator hits an extermal quark line it annihilates it, and 
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the inverse propagator with momentum of an adjoining internal quark line 
simply removes that propagator, replacing it by unity. The resulting "pinch" 
yields graphs such as Fig. WO 1 )- A pinch can change part of a vertex graph 
to a propagator graph, resulting in new terms that must be added to the 
conventional fermion propagator. Since the whole S-matrix element is gauge- 
invariant, it turns out that the new fermion propagator is gauge- invariant; 
gauge-dependent terms in the original Feynman graphs are cancelled by these 
other pinch terms. It is already known [5], [6J, [8 J that this procedure yields, in 
an i?£ gauge, a gauge-invariant gluon propagator of the type: 

KAk) = (<W " ^Hk) + ^ (62) 

where d(k) is completely gauge- invariant (independent of £) and the £ term 
receives no radiative corrections (other than wave-function renormalization). 
Not only do the longitudinal terms ~ k^ cancel out of the S-matrix, they 
also cancel out of the fermion pinch-technique propagator. (Because of this 
cancellation proper vertices obey naive ghost-free Ward identities.) Although 
I have illustrated the pinch technique only to one loop in the figure, it is pos- 
sible to extend it to all orders and to non-perturbative phenomena (see [16] 
and references therein). For the gap equation the essential non-perturbative 
phenomenon is a gluon mass, and I use d(k) = l/(k 2 + m 2 ). 

As is well-known [8], the pinch technique leads to Schwinger-Dyson equa- 
tions identical to those of the background-field method in the Feynman gauge 
(£ = 1). However, one can formulate the pinch technique in any gauge [UJ 
, in the sense that ghost-free Ward identities and certain other structural 
elements of Schwinger-Dyson equations important for the pinch technique 
in the Feynman gauge still hold. I will, in Pilaftsis' [47J sense of the pinch 
technique, use Landau gauge as all other workers do, because ultraviolet cor- 
rections to the vertex are unimportant to one loop in this gauge. Because 
the Wasd identity relating the vertex to the fermion propagator has no ghost 
terms, the same is true (as in QED) for the coefficient A(p 2 ) of pj in the 
fermion propagator [see Eq. (1301) ]. In the infrared regime A(p 2 ) should, like 
the running charge, not run much. So I expect A(p 2 ) ~ ^4(0) over a large 
momentum range, and will set it to unity. (With the gauge technique vertex, 
any constant A{p 2 ) cancels out from the gap equation.) 

The gauge-technique solution for a pinch-technique form factor is an 
infrared-valid approximation to the form factor that is asymptotically ex- 
act as the gluon momentum vanishes, and its ultraviolet inaccuracies can 



21 



be compensated for [IS], although I will not attempt that here. Because its 
Ward identity is ghost-free, just as in QED, all necessary formulas can be 
read off with straightforward modifications from the QED work in [18] . One 
pinch technique/gauge technique improper form factor is (omitting irrelevant 
group factors and indices): 

f,(p,p') = s(p)r„(p,p')s( P r ) = - r L^[( I j lll+1 ^)s(p>)-s(p)( lj i 1 ,+ 1 j)} 

(63) 

and obeys the QED-like Ward identity 

(p-p'),F ll (p,p') = S(p')-S(p). (64) 

Eq. (163]) is not a unique choice for the form factor, but it has the advantage of 
being an identity for a free massive theory, with T^p, p') = 7^ and S^ 1 (p) = 
p 1 — iM. All other choices for the pinch-technique form factor, as well as 
corrections needed for the ultraviolet behavior at loop level, are identically 
conserved and therefore differ from Eq. ( 1631) only by terms that vanish with 
at least one more power of p — p' at small values of p — p'. 
The Schwinger-Dyson (SD) equation is: 

S \P)S(P) = 1 + / d'k^k^p - k)F^p) (65) 

and with the gauge-technique form factor it is linear in the fermion prop- 
agator. In the pinch technique the gauge-boson propagator, which shows 
dynamical mass generation, has the form given in Eq. fl62l) . except that the 
gauge-dependent term is dropped. 

To study CSB I take S~ l {p) = p 1 and extract the coefficients of p 1 in the 
SD equation; this yields the gap equation. Because the form factor of Eq. (1631) 
is identically equal to that of a free massive field theory and because the 
running fermion mass is not running very fast in the infrared, the resulting 
pinch-technique equation at small momentum (including small integration 
momentum k) is really the same as I started with in Eq. (1291) . except that 
for small momenta I ignore the running of masses and the coupling. It is 
plausible that a nearly-correct modification of the infrared pinch-technique 
equation that accounts for ultraviolet effects is simply to let the masses and 
coupling run — cases I have already reviewed or analyzed. Actually, it is 
somewhat more complicated than this, but the final analysis is much too 
elaborate to discuss here. 
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The gauge technique incorporates infrared-important vertex corrections 
and the pinch technique yields a unique gauge-invariant gap equation in an 
arbitrary gauge for the underlying Feynman graphs. At the present rather 
simple level of approximation to the pinch technique/gauge technique gap 
equation, the deep-infrared gap equation, with plausible ultraviolet correc- 
tions, is the same as found by many others over decades by ignoring vertex 
corrections and working in the Landau gauge (of course, self-consistent at 
one-loop level). 



6 Finite temperature CSB 

The full finite-temperature (T) gap equation, especially for the pinch tech- 
nique, has a number of complications that I have not dealt with as of this 
writing, but will save for a more detailed work. These include separating the 
gluon propagator into space-space, space-time, and time-time components, 
which are no longer related by Lorentz invariance and which have two mass 
scales, the magnetic and electric (Debye). This complicates the application 
of pinch-technique cancellation mechanisms [12J. I will only present here 
the crudest initial attempts at understanding finite-temperature CSB with 
massive-gluon exchange, using what I call the superconductor approxima- 
tion, because it is in the same spirit as used in the original BCS paper on 
superconductivity. It amounts to saying that the gluon propagator is rela- 
tively unchanged by temperature effects as long as the temperature is not 
too large compared to the phase-transition temperature T c ~ 170 MeV. The 
reason is that the gluon already has a large T = mass of some 600 MeV, 
which is not changed drastically by thermal effects at T ~ T c . However, 
the fermion mass steadily decreases from its T = value, eventually van- 
ishing at the CSB phase transition^). The difference between fermionic and 
gluonic dynamical mass generation at finite T is that the effective coupling 
strength for fermions is decreased as T increases, by a factor something like 
tanh(/3a>/2) where (3 — 1/T and to is a characteristic energy scale; ultimately, 
the fermionic mass has to vanish as the coupling diminishes. But the gluonic 
mass is increased by a factor roughly coth(/3o;/2) (so the mass grows like T 
at large T). 

In the superconductor approximation I convert the original zero-temperature 

9 By mass I mean a quantity that violates chiral symmetry, which the usual perturbative 
finite-temperature fermion "mass" does not. 
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gap equation, Eq. (|2T)j) . to finite temperature by the usual replacement of the 
integral over the (Euclidean) time momentum fc 4 by a discrete frequency sum: 



/fr- £ ^4-2^ + 1)]. (66) 

N=— oo 



As mentioned above it is not quite correct simply to make this substitution 
in an equation such as (|2"5|) where the vector-propagator kinematics have 
been worked out at zero temperature. A further approximation is to average 
certain momentum- dependent quantities that vary fairly slowly by averages, 
which allows us to make contact with an already-studied zero-T equation, 
Eq. ( l44l) . The resulting gap equation is then Eq. ( 1441) with the modified 
zero-momentum but finite-T coupling: 

g \k = 0,T = 0) - g\k = 0,T)(tanh(^)> = G 2 (T) (67) 

with the finite-T coupling determined from the zero-momentum form of 
Eq. (I40p by using a plausible form for the temperature dependence of the 
finite-T gluon mass: 

g\k = 0,T) = {b\n[^^l]}- 1 ; m(T) = m(T = 0) coth(^) (68) 

where I choose for the average gluon frequency 

47rm(0) 

UG = iV(fc = 0,T = 0) (69) 

so as to give a simple but fairly accurate high-T limit ( cf the value used 
in connection with the FSE of m(T) = 0.16Ng%, with g\ ~ g 2 (0)T). The 
fermion frequency uop is of the form: 



uj F = y/k 2 + M 2 (k). (70) 

At the level of the simple approximations used here, it is not possible to 
predict with any accuracy the actual ratio of the temperature T x a, at which 
adjoint chiral symmetry breaking is restored, to the usual deconfinement 
transition temperature T^. However, after some not very interesting numer- 
ics involving plausible ranges of unknown and approximated quantities, and 
using the superconductor approximation given here, it appears possible that 



24 



T x a > T d . If this possibility survives more detailed numerics then adjoint 
chiral symmetry breaking survives deconfinement, as the lattice simulations 
show. 

Of course, there is much more to learn about QCD in general from adjoint 
CSB at finite temperature beyond our interests here, as the literature shows 
[231 [2U ES] ■ I hope to take up these questions in more detail later. 

7 Summary 

I use the FSE to estimate, using d = 3 dynamical mass calculations, the 
usual d — 4 strong coupling at zero momentum: a s (0) ~ 0.5 (for three light 
flavors), a value in accord with other estimates using the pinch-technique 
Schwinger-Dyson equations. This is somewhat smaller than phenomeno- 
logical estimates of around 0.75. I then consider modern versions of the 
fermion gap equation both for quarks and for adjoint fermions in QCD— 
made modern by adding dynamical gluon mass effects and by consideration 
of the gauge-invariant pinch techniqe/gap technique gap equation. The final 
results, although not impressively accurate, suggest that confinement is es- 
sential to break CSB for quarks, but that standard fermion gap physics can 
explain CSB for adjoint fermions, for some (uncalculated) range of tempera- 
tures, possibly reaching above the deconfinement temperature. These results 
are consistent with present-day lattice simulations. 
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